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ABSTRACT 

Recent years have shown many advances in our knowledge of the collisional evolution of proto- 
planetary dust. Based on a variety of dust-collision experiments in the laboratory, our view of the 
growth of dust aggregates in protoplanetary disks is now supported by a deeper understanding 
of the physics involved in the interaction between dust agglomerates. However, the parameter 
space, which determines the cohisional outcome, is huge and sometimes inaccessible to labora- 
tory experiments. Very large or fluffy dust aggregates and extremely low collision velocities are 
beyond the boundary of today's laboratories. It is therefore desirable to augment our empir- 
ical knowledge of dust-collision physics with a numerical method to treat arbitrary aggregate 
sizes, porosities and collision velocities. In this article, we implement experimentally-determined 
material parameters of highly porous dust aggregates into a Smooth Particle Hydrodynamics 
(SPH) code, in particular an omnidirectional compressive-strength and a tensile-strength rela- 
tion. We also give a prescription of calibrating the SPH code with compression and low-velocity 
impact experiments. In the process of calibration, we developed a dynamic compressive-strength 
relation and estimated a relation for the shear strength. Finally, we defined and performed a 
series of benchmark tests and found the agreement between experimental results and numerical 
simulations to be very satisfactory. SPH codes have been used in the past to study collisions at 
rather high velocities. At the end of this work, we show examples of future applications in the 
low-velocity regime of collisional evolution. 

Subject headings: accretion, accretion disks - methods: laboratory, numerical - planetary systems: for- 
mation - solar system: formation 

1. Introduction 

1.1. Protoplanetary dust growth 

The formation of planetesimals, the km-sized 
solid bodies whose further growth is controlled by 
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mutual gravitational attraction, is still enigmatic. 
Collisions among the dust aggregates are con- 
trolled by Brownian motion, drift motions with 
respect to the gas of th e protoplanetary disk , 
and turbulence in the gas ( Weidenschilling 19 77al : 
Weidenschilling and Cuzzil 19931 ). Once in con- 
tact, two dust gr ains experience a mutual van 



der Waals force ( Heim et aD Il999[ ). From the 



theoretical and experimental standpoints, it is ev- 
ident that the (sub-)micrometer-sized protoplan- 
etary dust grains initially undergo hit-and-stick 
collisions, whic h lead to the formation o f frac- 
tal aggregates (IWeidenschilling and Cuzzil I1993I : 

l2004h . As 



Blum et all I2QQQI : iKrause and BlumI l2QQ4f ). As 
the collision energy increases, due to increasing 
aggregate mass and collision velocity, dust aggre- 
gates undergo a restructuring phase, i n which they 
acqu ire denser structures (IPomini k and TielensI 



1997; 
2008; 
iments 



Weidling et all I2009D . 



Blum and WurmI I2OOOI : IWada et all 12007 



Laboratory exper- 
showed that collisions among the dust 
aggregates result in fragmentation, i.e. in mass 
loss, if the impac t velocities exceed ~ lms~^ 
( Blum and WurmI [2OO8). Depending on the 
disk model, this means that the direct colli- 
sional growth process ends (at the latest) at 
aggregate sizes for which this velocity is ex- 
ceeded. For a minimum-mass solar nebula model 



(|Weidenschillindll977bl : iHavashi et al.lll985h . this 
size is approximately 10 cm. 

The further gro wth is still highl y speculative. 
Wurm et all (|200lh and lBlumI (|2004h proposed the 
accretion of collisional fragments by aerodynamic 
and el ectros tatic effects, respec t ively. I Wurm et al. 



(|2005h and iTeiser and WurmI (|2009h showed ex- 
perimentally that a fraction of a dust projectile 
can stick to a solidified larger dust target even 
at very large velocities. None of these processes, 
however, seem to work globally and under all cir- 
cumstances so that very specific conditions are re- 
quired for the dust aggregates to grow at high im- 
pact velocities. There is clearly a lack of under- 
standing the detailed physics involved in the colli- 
sions between macroscopic dust aggregates of ar- 
bitrary composition and porosity. Without better 
knowledge of the collisional physics of these bod- 
ies, any attempt to model the formation of plan- 
etesimals as an aggregation process will have to 
fail. 



1.2. Previous work 

In the three previous papers of this series, we 
described the collisional physics of high-porosity 
protoplanetary dust aggr egates up to the cm- 
size regime. In paper I ( Blum et al.l I2OO6 ). we 
introduced a method to experimentally produce 
monolithic dust aggregates with diameters of 2.5 
cm. By choosing either monodisperse spherical 
monomer particles, quasi-monodisperse irregular 
particles, or polydisperse irregular grains, we pro- 
duced dust aggregates with volume filling factors 
(i.e. packing densities) (j) = p/ ps of (j) = 0.15, 

= 0.11, and (f) = 0.07, respectively (see Table 

1 in paper I for more details about the monomer- 
particle properties). Here, p and ps are the ag- 
gregate and the monomer density. Static uniaxial 
compression of these dust samples revealed that 
the maximum compaction for these high-porosity 
dust aggregates is (/)max = 0.20 ... 0.33, a value 
very close to the overall porosity found in comets. 
The tensile strengths of our dust samples were de- 
termined to |T| = 200 ... 6, 300 Pa, depending on 
the monomer properties and the compaction. Also 
these val ues are close to those f ound for comets. 
Paper II ([Langkowski et al.ll2008h concentrated on 
low-velocity impacts into these high-porosity dust 
samples. We showed that sticking by penetra- 
tion is the dominating process for impacts above 
a threshold velocity of ~ 1 ms~^ for projectiles in 
the mm-size regime and flat dust targets. For shal- 
low penetration, i.e. for impacts below the thresh- 
old velocity, the projectiles bounce off, leaving a 
well-defined crater. It is obvious that the colli- 
sions resul t in the compaction of the target. In 
paper III dWeidling et a ll 120091) . we investigated 
the compaction for high-porosity mm-sized dust 
aggregates in bouncing collisions. Bouncing colli- 
sions amon g dust aggregates show c onsiderable en- 
ergy losses ( Blum and Miinchll 19931 ) so that it was 
natural to assume some degree of compaction. In 
paper HI, we found that - although a single colli- 
sion leads only to very localized compaction of the 
dust aggregate - mm-sized dust aggregates in pro- 
toplanetary disks can reach volume filling factors 
of ^ 0.35 within a few dozen years. 

1.3. Objectives 

All previous experiments (see Sects. 11.11 and 
II. 2p were limited by the experimentally avail- 
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able dust- aggregate sizes and morphologies and 
the achievable collision velocities. In the astro- 
physical context, the need for numerical simula- 
tions of collisions between dust aggregates of arbi- 
trary composition, size and impact velocity arises 
from the fact that only a limited parameter space 
can be covered by experiments. The ongoing de- 
bate about threshold velocities for sticking, bounc- 
ing, compaction, and fragmentation as well as the 
fragment size distribution requires a thorough in- 
vestigation of a wide range of collisions, varied 
over supposedly critical parameters, such as col- 
lision velocity, porosity, size, impact parameter, 
impact angle and shape of the colliding dust ag- 
gregates. An extensive parameter study of that 
kind is not feasible under laboratory conditions 
for the parameter ranges in question. Therefore, 
we aim to calibrate a Smooth Particle Hydrody- 
namics (SPH) code and validate this model thor- 
oughly with a series of independent benchmark 
tests. Hence, the SPH code gains a deeper re- 
liability and the conducted numerical simulations 
provide well-grounded insight into the physical be- 
havior of dust aggregates. 

2. SPH in Dust Collisions 

SPH is a meshless Lagrangian particle method 
originally developed for astrophysical hydrody- 
namics applications. A detailed description of 
the original SPH method may, e.g., be found in 
[Monaghan (2005). The SPH code we utilize for the 
simulations in this work and the underlying poros- 
ity model are intro duced and de scribed in full 
depth in [Ceretshaus er et al.l ( 20091 ). In the 1990s, 
SPH has been extended to model the elastic and 



plastic beha vior of solids, see e.g. iLiberskv et al 



(|1993D and iRandles and Liberskvl (|1996D . The 
continuous solid objects are discretized into inter- 
acting mass packages called particles, which form 
a natural frame of reference for any deformation 
and fragmentation that may occur. 

The SPH code solves the equations of contin- 
uum mechanics in Lagrangian form, in particular 
the continuity equation 



dt ^ ^ dxa ' 
and the equation of motion 

dVg _ 1 daaf3 

dt p dxp 



(1) 



(2) 



Here, Einstein's summing convention holds through- 
out and Greek indices denote spatial coordinates. 
The variables have their usual meanings, i.e., p 
denotes the density, v the velocity, and aaf3 the 
stress tensor. The latter is defined according to 



(3) 



consisting of a pressure part with pressure p and 
and a shear part given by the traceless deviatoric 
stress tensor Sap- 

The deviatoric stress is defined by the consti- 
tutive equations. To model elastic behavior ac- 
cording to Hooke's l aw we adopt the approach by 
Benz and Asphau^ (1994) for the time evolution 
of the deviatoric stress, 

(4) 

where /i is the shear modulus and d denotes the 
dimension. The rotation rate tensor Raf3 reads 



Rex (3 



1 / dva_ _ dvp\ 



2 \dXf3 dxaj 
and the strain rate tensor €^[3 accordingly 



2 \dxp dxa J ' 



(5) 



(6) 



This set of equations is closed by a suitable equa- 
tion of state and describes the elastic behavior of 
a solid body. 

In order to simulate also the plastic behavior of 
porous bodies, we adopt a mod ified version of the 
porosity model bv lSironol (|2004f ) (Fig.[l]). Accord- 
ing to this approach, plasticity is modeled within 
the equation of state, which is divided into three 
different regimes. In the first regime, plastic be- 
havior is caused by compression that exceeds a 
critical limit, the compressive strength S(p), while 
in the second regime, tension exceeds the ten- 
sile strength limit T(p). In between these lim- 
its, the third, the elastic regime of the material is 
described by a special version of the Murnaghan 
equation of state. Thus, the full equation of state 
reads 



P{p) 



S(P) 

K{p'o){p/p'o 
T{p) 



1) Pc 



P> Pt 
' <P<_pt 

P< Pc 



(7) 
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Fig. 1. — The modified Sirono porosity model is 
divided into the elastic, compressive and tensile 
regimes. The latter two are reached by exceeding 
the compressive and tensile strength, respectively, 
which leads to irreversible plastic deformation. 

The quantity p' ^ denotes the reference density, 
which is the density of the material without any 
external stress. p+ and p~ are limiting quanti- 
ties, where the transition between the elastic and 
plastic regime for compression and tension, respec- 
tively, takes place. Once these limits are exceeded, 
the material leaves the elastic path which repre- 
sents the path where energy is conserved, and loses 
internal energy by following the paths of the com- 
pressive and tensile strength (Fig. [T]). 

3. Towards an Equation of State for Dust 
Aggregates 

In this laboratory section, we will provide the 
macroscopic material parameters, which are nec- 
essary for the SPH model introduced in section [2l 
We recapitulate the tensile str ength measurements 
of lBlum and Schraplerl ( 20041 ) and give an interpo- 
lation for different volume filling factors. The com- 
pressive strength for unidirec tional (ID) compres- 
sion w as also measured by um and Schraplerl 
(|2004t ). while in this paper we will present mea- 
surements on omnidirectional (3D) compression. 
Moreover, we will introduce a simple impact ex- 
periment, which will be used for calibrating the 
SPH model: A mm-sized glass bead (or a glass 
bead analog) impacts into a well-defined 2.5 cm 
dust sample at a collision speed of 0.1 to 1 m s~^. 
The dust sample, consisting of 1.5 /im Si02 



monodisperse spheres, was formed by random bal- 
listic deposition (RBD) and h as therefore a volume 
filling factor of 0o = 0-15 (see lBlum and Schrapler 
(|2004l ) and references therein). The deceleration 
curve, penetration depth and impact duration of 
the glass bead are measured as well as the com- 
pression of the dust beneath the glass bead to 
compare these results with an impact computed 
by the SPH model. 

3.1. Tensile Strength 

In Blum and Schraplerl (2004) and paper I, we 



reported on measurements of the tensile strength 
of dust samples of various constitutions (i.e. 
monomer size distribution, morphology, compo- 
sition, volume filling factor). The best set of data 
was collected for the dust aggregates consisting 
of spherical 1.5 /im Si02 monomers (see above). 
For packing densities of = 0.15, = 0.41, 
(\) = 0.54 and = 0.66, we found tensile strengths 
of |T| = 1,000 Pa, |T| = 2,400 Pa, |T| = 3,700 
Pa and |T| = 6,300 Pa, respectively. To a good 
approximation, these values can be expressed by 
a relation of the form 



= - {m 



,2.8+1.48 (; 



(8) 



This expression will be used throughout this paper 
for the packing-density-dependence of the tensile 
strength. 

3.2. Static Measurement of Compressive 
Strength Curves 

The compression curve of a given material tells 
us how the material behaves under an applied 
pressure S in changing its volume filling factor cj). 
If the material can be described by macroscopic 
parameters, the volume filling factor is represen- 
tative for the material density and so the devel- 
opment of the compression curve ^(S) (cf. Eq. [7j) 
is essential to establish a collision model and learn 
about collisions of protoplanetary dust aggregates. 

Measurements of t he compression curve were 
already performed bv iBlum and Schraplerl ( 20041 ) 
and in paper I. We will again focus on the 
dust samples made of 1.5 pm Si02 spheres, 
whose properties a r e com piled in Table 1 in 
Blum and Schrapler (l2004h. In th e coni pression 
experiments of Blum and Schraplerl (2004), a dust 
sample was fixed between two parallel glass plates. 
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which were then pushed together with an increas- 
ing force. The measurement of the dust mass, 
dust volume, compression force, and, thus, pres- 
sure yield the compression curve ^(S). The force 
was applied in one direction, which is therefore 
called unidirectional compression. The dust sam- 
ple flattens in the direction of the force but, at the 
same time, also expands in the other directions. 
For dust samples made of 1.5 fim Si02 spheres, 
this leads to an equilibrium fllling factor of 0.33 
for pressures exceeding 10^ Pa. This compression 
curve is only applicable to protoplanetary dust 
collisions, if the material compressed in the im- 
pact zone creeps sideways as it did in the static 
experiments. As we will show later by x-ray anal- 
ysis of the compression next to an impact site, 
this is not the case. 

Consequently, a second way to measure the 
compressive strength curve is to flx the dust sam- 
ple at the sides with closed walls. In this case, the 
pressure cannot be released and acts from all sides, 
thus omnidirectional compression. We performed 
experiments in which we cut a cylindrical section 
from an RBD dust sample with a thin- walled plas- 
tic tube of 7 mm diameter. This cylindrical dust 
sample of approximately 1 cm height was then put 
into a 7 mm borehole in an aluminum block. Care- 
fully pushing a piston into this borehole leads to 
an omnidirectional pressure onto the dust sample 
(see inset in Fig. [2]). The setup was put onto a 
balance and the piston was loaded with weights of 
increasing mass. This weight force, divided by the 
piston area, yields the pressure E, while the mass 
and height of the dust sample determine the vol- 
ume fllling factor (j) (Fig. [2j). Due to the fact that 
the dust sample is not a frictionless fluid, force 
chains inside the sample might locally reduce the 
pressure. Thus, the pressure for the idealized com- 
pression curve can be slightly lower. 

The solid line in Fig. [2] denotes an analytical 
approximation of the mean fllling factor of nine 
individual experiments as a function of the ap- 
plied pressure and the gray shaded area is the 
standard deviation of the measurements. The an- 
alytical function is based on a Fermi distribution 
with logarithmic pressure in the energy term 



(9) 



exp (^k^^gPHLj + 1 
and is only valid for (j) > (po- Foi" pressures below 



Table 1: Parameters of the analytical approxima- 
tion for the two compression curves 

^2 Pm [kPa] A [dex] 

unidirectional 0.15 0.33 5.6 0.33 
omnidirectional 0.12 0.58 13.0 0.58 



11(^0 ) the dust aggregate behaves elastically. The 
parameters for the unidirectional and omnidirec- 
tional compression curve are given in Table [H The 
bottom plot in Fig. [2] gives the deviation between 
the analytical approximation and the data, which 
is within 0err = ±0.02. Often, the inverse function 
i;(^) is used (see Eq. [7]), which is here 



S((/)) = 



1 



A-lnlO 



(10) 



Compared with the uni directional compres- 



sion c urve (dashed line) of iBlum and Schrapler 
(|2004l ). the fllling factor also starts off at the orig- 
inal dust sample fllling f actor of (j)o = 0.15 (cf. 
Blum and Schraplerll2004l ). but diverges from the 



unidirectional curve for pressures p ^ 10 Pa. 




1 0^ 1 04 1 0^ 
pressure [Pa] 



Fig. 2. — The omnidirectional (solid line) and uni- 
directional (dashed line) static compression curves 
and the standard deviations of the measurements 
(gray shaded areas). The upper plot shows the an- 
alytical approximations from Eq. [9] and Table [H 
while the lower plot gives the deviation between 
approximation and measurement. The inset de- 
scribes the setup for the omnidirectional measure- 
ment. 
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For those pressures, the fihing factor is systemat- 
ically higher, raeaning on the other hand that the 
same filling factor is much easier to achieve if the 
pressure acts from all sides. So far there was no 
equilibrium filling factor found like in the case of 
the unidirectional compression experiments. The 
filling factor still significantly increases for the 
highest applied pressure of 10^ Pa. However, the 
analytical approximation indicates an equilibrium 
for 02 = 0.58, which is not far from random close 
pac king of monodisperse spheres (0 c:^ 0.64, see 
e.g. IXorquato et al.l HoOOh . the maximal possible 



release 
mechanisiTi 



compression without breaking the dust grains. 

The new compression curve is still a static mea- 
surement. It is applicable for omnidirectional 
static pressures like the hydrostatic equilibrium 
inside planetesimal bodies. It is questionable if 
this compression curve is valid for dynamic colli- 
sions but it is a second attempt to assume that 
surrounding material, which does not interact in 
a collision, acts as a confining wall to the active 
impact volume instead of creeping sideways. 

3.3. Deceleration Experiments 

3.3.1. Experimental Setup 

The experimental setup consists of a vacuum 
chamber (gas pressure ~ 0.5 mbar) in which a pro- 
jectile is suspended on a thin fiber (Fig. [3|) with 
negligible mass to prevent rotation and lateral ve- 
locities. The distance between the suspended pro- 
jectile and the surface of the dust sample deter- 
mines the impact velocity vq- The projectile con- 
sists of an elongated solidified epoxy droplet at 
the bottom and a cylindrical plastic tube at the 
top end. After the release of the projectile, it is 
accelerated by gravity and decelerated once it is 
in contact with the dust sample. The deceleration 
within the dust sample is observed by a high-speed 
camera (Fig. H]). From the deceleration curve of 
the projectile we can derive fundamental dynamic 
properties of the target dust aggregate. 

The bottom shapes of the projectiles were 
spherical with diameters of R ~ 0.5 mm and R ^ 
1.5 mm and masses of m ~ 1 mg and m ~ 30 mg, 
respectively (see Table [2]). The effective densities 
of the projectiles of p = 2, 400 . . . 3, 100 kg m-^ 
match those of the astronomically relevant sili- 
cates, while the combination of low-density epoxy 
and plastic tube made sure that the top of the 



fiber 



high-speed 
camera 



i 



projectile 
(glass bead 
analog) 




dust sample 



Fig. 3. — Setup for the deceleration measurement: 
An elongated projectile as a glass bead analog was 
dropped into the dust sample. Before dropping 
from a height of 1 to 40 mm, it was suspended 
on a fiber with negligible mass to avoid rotational 
motion. A high-speed camera observes the decel- 
eration of the projectile. 




Fig. 4. — Image sequence of a decelerated projec- 
tile. The time between two images is 1.7 ms. 

projectile was always visible to the camera even if 
the intrusion depth was larger than the projectile 
diameter. The high-speed camera was operated 
at a frame rate of 12,000 frames per second with 
a resolution of ~ 30 /im/pixel and the position 
of the upper edge of the projectile was measured 
with sub-pixel accuracy of ~ 3 /im. The first 
touch of the projectile with the surface of the dust 
sample marks the time t = and can clearly be 
determined from the deviation of the trajectory 
compared to a free falling projectile. After its 
deepest penetration, the projectile bounces back 
(by ^ 100 /im) and oscillates in the vertical direc- 
tion, which we will not take into account in the 
further discussion. 
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Table 2: Experimental results of the deceleration experiments. 


experi- 


projectile 


projectile 


effective 


impact 


penetration 


stopping 


standard 


standard 


ment 


diameter 


mass 


projectile 


velocity 


depth 


time 


deviation a for 


deviation a for 


number 


2 • R [mm] 


m [mg] 


density [kg/m*^] 


r — It 


D [mm] 


T [ms] 


polynomial [/im] 


sine [/im] 


1 


2.73 


25.7 


2412 


0.89 


3.16 


5.92 


2.43 


21.48 


2 


2.94 


32.0 


2404 


0.85 


3.08 


5.93 


2.99 


22.64 


3 


2.77 


26.4 


2372 


0.73 


2.70 


6.09 


1.87 


17.29 


4 


2.94 


32.0 


2404 


0.17 


0.70 


6.18 


4.63 


10.78 


5 


2.94 


32.0 


2404 


0.16 


0.75 


7.13 


2.84 


7.69 


6 


2.94 


32.0 


2404 


0.20 


0.80 


5.92 


2.96 


8.95 


7 


2.94 


32.0 


2404 


0.32 


1.31 


6.47 


1.62 


2.53 


8 


2.77 


27.6 


2480 


0.50 


2.03 


6.38 


1.82 


3.57 


9 


2.77 


26.4 


2372 


0.37 


1.33 


5.71 


3.65 


6.26 


10 


0.99 


1.5 


2854 


0.67 


1.06 


3.09 


3.48 


27.77 


11 


0.99 


1.5 


2854 


0.76 


1.15 


2.94 


1.41 


23.24 


12 


0.99 


1.5 


2854 


0.79 


1.47 


3.12 


14.29 


5.62 


13 


0.85 


1.0 


3109 


0.19 


0.32 


2.46 


1.94 


3.79 


14 


0.85 


1.0 


3109 


0.35 


0.83 


3.28 


2.60 


11.15 


15 


0.85 


1.0 


3109 


0.36 


0.72 


3.11 


2.71 


3.73 



3.3.2. Experimental Results 

We performed 15 impacts of our projectiles 
into the porous dust samples, which are compiled 
in Table [2l The time-resolved deceleration data 
h{t) were cleaned from gravitational influence by 
adding ^gt'^ to the negative intrusions so that the 
gravity-independent deepest penetration depth D 
and stopping time T could be determined. The in- 
trusion curves were normalized in space and time 
through h\f) = h{t)/D and t' = t/T so that 
h\t' = 0) = (first contact) and h\t' = 1) = -1 
(deepest intrusion), and can then be well repre- 
sented by a sine function 

/i'(t') = -sin(t'.|) . (11) 

Alternatively, a fourth order polynomial with 
only one free parameter was used for fitting 
the data, where the mean standard deviation 
between the fit and the TV data points a = 

\l W Xlili(^'(^i) ~ ^i)^ amounts to only 2-4 /im 
in absolute units (cf. Table [2]). Although the 
standard deviation for the sine function is rather 
of the order of 10 /im, we take the sine function 
because it has no free parameter and the standard 
deviation is still less than the pixel size of 30 /im. 
In few experiments, the 1 mm projectiles canted 
over before coming to rest. In these cases, the 
data was used as long as reliable and the remain- 
ing deceleration curve and, thus, the penetration 
depth was extrapolated. 

Figure [5] shows all measured deceleration curves 
in absolute units. Different intrusion depths and 
stopping times can clearly be distinguished in this 



plot. The intrusion depths increase with increas- 
ing impact velocities (i.e. with the absolute values 
of the initial slopes of the curves), while the stop- 
ping times are rather constant for one projectile 
size (T ~ 6 ms for 3 mm projectiles [nos. 1-9 in 
Fig. [5] and T ~ 3 ms for 1 mm projectiles [nos. 10- 
15 in Fig. [5]) and, thus, independent from the im- 
pact velocity vq. A test yielded the best-fitting 
power-law relation between the penetration depth, 
impact velocity and mass of the form 

D = jB- • , (12) 

with QfD = 0.23±0.13, (3b = 0.89±0.34, and 7d = 
(3.86dz0.11)-10-2kg-0-23m0-iis-0-^^ (Fig.[6]). The 
respective errors denote the la uncertainties. A 
more intuitive relation would he D (x mvoA~^^ 




2 4 6 

time [ms] 



Fig. 5. — Fitted deceleration curves in absolute 
units. The experiment numbers at the decelera- 
tion curves confer to those in Table O 
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0.01 0.10 

^0.23,, 0.89 ri>^0-23/^ /„\ 0.891 

m Vo [kg vrn/sj J 

Fig. 6. — The best relation for the penetration 
depths from a test yields a dependence of I) ex 
^0.23^0.89 rpj^g intuitive relation D oc mvA~^ is 
possible within the uncertainties. 




Fig. 7. — Dust sticks to the projectiles after 
pulling them out of the dust sample. This is an 
indication of compacted material under the pro- 
jectile as will be confirmed in Sect. 13. 4[ 

with A = 7tR^ being the cross section of the pro- 
jectile. This relation has a clear physical mean- 
ing as the penetration depth is determined by the 
quotient of the momentum mvo as driving force 
and the cross sectional area A as resistive pa- 
rameter. With ao = I and = 1 being pos- 
sible within the uncertainties, the linear relation 
D oc mvoA~^ oc pp^^Rvo is also possible. How- 
ever, constraining the exponent as I) oc p^^ 
in Eq. [12] was unfortunately not feasible due to 
the too small variations in the effective projectile 
density pp (cf. Table [2]). 

For the stopping time, we found 

T = 7t • m"^ • v^^ (13) 

with ar = 0.23 ± 0.08, f3T = 0.01 ± 0.23, and 
7T = (6.77 ± 0.20) • 10-2 ]^g-o.23^-om^im 



found for the stopping time. 

After pulling the projectiles out of the dust 
sample, dust stuck to the surface with which it 
had been in contact before (Fig.[7j). With the pre- 
liminary assumption that this is compacted dust 
and the layer where it broke off is the transition 
from compacted to non-compacted dust (transi- 
tion in tensile-strength), this gives an indication 
for the compressed volume which will be analyzed 
in detail in the forthcoming section. 

3.4. Dynamic Compression Experiments 

3.4' Experimental Setup 

In order to investigate in more detail the com- 
pression behavior of the dust aggregates by colli- 
sions, we performed impact experiments with sub- 
sequent x-ray micro-tomography measurements to 
analyze the degree of compaction. 

Under vacuum conditions we dropped a single 
glass spherule with a diameter of ~ 1 mm from a 
given height of ~ 75 mm into an RBD dust sample 
within a plastic tube with 7 mm diameter. To en- 
sure the sphere to preferably hit the center of the 
dust sample within the narrow plastic tube, the 
released projectile was guided by falling through a 
tube. Due to friction and collisions with the tube's 
walls the impact velocity of (0.8±0.1) m s"^, that 
was independently measured by high-speed imag- 
ing in 10 drops, is much lower than expected from 
free fall. However, the velocity in the specific ex- 
periment was not measured and can well be in 
the lower range of the error. From comparison of 

X-ray source detector 




sample carrier 

Fig. 8. — Setup of the x-ray micro-CT measure- 
ment: the sample is rotated between an x-ray 
source and a detector. A 3D density reconstruc- 
tion can be computed from the transmission im- 
ages. 
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the observed penetration depth (see Fig.[9]in Sect. 
I3.4.2P with the results in Fig. [6] we expect a veloc- 
ity of V = 0.65 m s~^, which we will use for the 
further study. 

For analyzing the density distribution of the 
dust sample cutout with the embedded glass 
sphere, the dust sample was scanned by an x-ray 
micro-computer-tomograph (Micro-CT SkyScan 
1074) at the University of Osnabriick. The dust 
sample was positioned on top of a rotatable sam- 
ple carrier between the x-ray source and the de- 
tector (CCD camera) (Fig. [8j). While rotating 
stepwise around by 360°, 400 transmission im- 
ages were captured. Based on this data set, a 
3-dimensional density reconstruction was calcu- 
lated by the SkyScan Cone-Beam Reconstruction 
Software provided with the x-ray micro- CT in- 
strument. 

3.4-2. Experimental Results 

In the following we present the results of two 
impact experiments. Further experiments with 
differently sized spheres and different impact ve- 
locities are intended. To visualize the spatial den- 
sity distribution of the observed dust sample with 
the impacted glass sphere, the 3-dimensional re- 
construction data was cylinder-symmetrically av- 
eraged with the vertical axis aligned with the 
sphere center. Figure [9] displays the mean volume 
filling factor as a function of height and radius, 
whereas the data is mirrored with respect to the 
vertical center line of the diagram. The color gra- 
dient from yellow to light blue underneath the im- 
pacted sphere (red color: saturated density values 
of the considerably denser glass spherule) clearly 
shows the densification of the porous dust sam- 
ple with an initial volume filling factor of (j)o ~ 
0.15. The compressed area, emphasized by over- 
plotted contour-lines, is located almost cylindri- 
cally shaped beneath the sphere and extends only 
slightly to the lateral borders of the sphere. Thus, 
the assumption of an omnidirectional compression 
curve, made in Sect. 13. 2[ seems to be justified. 

Analysis of the distribution of occurring volume 
filling factors related to their fraction of volume 
within an uncompressed dust sample provides a 
Gaussian-shaped distribution with a mean value of 
(j) ^ 0.15 (Fig. [TO]). Figure HU (top) shows the vol- 
ume fraction (normalized by the sphere volume) of 
volume filling factors, which we determined only 
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Fig. 9. — Spatial averaged volume filling factors 
mirrored at the vertical center line. The volume 
under the sphere is compacted to a volume fill- 
ing factor > 0.2 (yellow), while the surrounding 
material is nearly unaffected (light blue). 
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Fig. 10. — Distribution of volume filling factors 
for an uncompressed dust sample, which follows a 
Gaussian distribution with a mean value of ~ 
0.15. 

regarding the compacted volume underneath the 
impacted sphere for the two impact experiments. 
In both curves the most prominent volume fill- 
ing factor is around (j) = 0.23, indicated by the 
dashed and dotted lines. The decreasing left flank 
of the curves corresponds to the transition region 
between compressed and uncompressed dust ma- 
terial (see right curve flank of Fig. [TO]). The same 
data plotted in a cumulative way (Fig. (TTJ bot- 
tom), represent the amount of compacted volume 
in units of the sphere volume that complies with a 
volume filling factor greater than a certain value. 
According to the volume filling factor values at 
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Fig. 11. — Top: Distribution of volume filling 
factors only for the compressed area underneath 
the impacted sphere for two experiments. The 
dashed and dotted lines mark the most occurring 
volume filling factors for each curve, lying at (j) ^ 
0.23. Bottom: Normalized volume fraction of 
compacted area corresponding to a volume filling 
factor > 6. 



the boundary to the uncompressed dust, given by 
the minima of the left side of the curves in Fig. [11] 
(top), we can conclude from the cumulative curves 
(Fig.[TTl bottom) that the compressed volume due 
to an impacting sphere of 1 mm size into a high- 
porosity dust sample {(j) ^ 0.15) fills the volume 
of ~ 0.8-1.2 sphere volumes. 

3.5. Requirements of a Dynamic Com- 
pressive Strength Curve 

As seen in the previous sections, we have abun- 
dant indirect information about the compression 
behavior of loose dust samples. However, the basic 
question how the dynamic compressive strength 
curve, ^(S), looks, remains unanswered. We ap- 



proach this problem the following way: (1) For 
low compressions, E ^ 0, the volume filling fac- 
tor is given by the initial properties of the mate- 
rial, i.e. (j) ^ (j)! (see Table [T]). (2) The maximum 
compression for S ^ oo is given by the value 02 
in Table [T] for the omnidirectional case, because 
the XRT analysis shows no material creeping side- 
ways a^_wa^_tlie_case_Jor fiow 
(see iBlum and Schraplerl 12004 ). (3) With these 
two limits in mind, we apply Eq. [9] as an ap- 
proximation to the functionality of the dynamic 
compressive strength, which leaves us with the 
two free parameters A and Pm- The maximum 
slope of the compression at E = is given by 
d^/dlgS = (02 — 0i)/A. For the unidirectional 
and omnidirectional static curves, we get slope val- 
ues of 0.55 and 0.79, respectively (see Table [T]). 
These are in fact not so different so that we adopt 
for the dynamic case the slope of the omnidirec- 
tional compression. Thus, we assume A = 0.58 
dex for the dynamic case. A refined study that 
takes both, A and as free paramet ers will 
be conducted in lGeretshauser et al.l ( 20091 ). but in 
this paper we only vary . 

4. Calibrating the SPH code 

The laboratory experiments in the previous sec- 
tion provided the static omnidirectional compres- 
sive strength E and the tensile strength relation 
T as most important ingredients for the Sirono 
porosity model implemen ted in the SPH code by 
iGeretshauser et aJ] (|2009l ). However, as it was al- 
ready pointed out in the laboratory section, the 
compressive strength relation has to be consid- 
ered dynamically. The only free parameter (see 
Sect. 13. 5p cannot be determined by experiments. 
Hence it has to be constrained by a numerical pa- 
rameter study. We will use the stopping time of 
the impacting glass bead as reference for this pa- 
rameter. 

In addition, a relation for the shear strength 
is also very hard to measure in the laboratory. 
Therefore, we suggest three simple relations de- 
pending on the dynamic compressive strength and 
tensile strength relations and use the qualitative 
comparison of the filling factor profile under the 
glass bead after impact to constrain this unknown 
quantity. 

Finally, we utilize the remaining experimentally 
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measured independent features of the experiments 
described in the laboratory section to vahdate our 
cahbration. 

4.1. Benchmark test - setup 

The given experimental setup (see Sect. 13. 4p 
was modeled with high resolution in two dimen- 
sions. Initially, the SPH particles were put on a 
triangular grid. All simulations were performed 
with the influence of gravity taken into account. 

The projectile was modelled with a circle of 
1.1 mm in diameter consisting of 1519 SPH parti- 
cles. Its material properties were simulated using 
the Murnaghan equation of state 



which is c = 30 m s ^ ( 


Blum and WurmI 20081: 


Paszun and DominikI 2008 


). 



p ■ 



Ko 

n 



^1 -1 

Po 



(14) 



with po = 2540 kg m ^ (total 2D mass per unit 
length m2D = 2.4-10-^ kg m"^), Kq = 5.0-10^ Pa 
and n = 4. The density has been chosen such that 
it matches the experimental specifications. The 
other material parameters are similar to those of 
sandstone. They can be found to gether with th e 
Murnaghan equation of state in MeloshI (|l989l ). 
The exact choice of the bulk modulus Kq and the 
Murnaghan exponent n does not have significant 
effects. The glass bead was treated as fully elastic. 
The impact velocity was 0.65 m s~^. 

The dust sample was modelled as a 8 x 5 mm^ 
rectangle with 64421 SPH particles. About 
0.15 mm at the bottom and 0.56 mm at each side 
of the rectangle were used as refiecting bound- 
ary by setting their acceleration to zero at each 
time step. The porous material was simulated by 
using the modified version of the Sirono model 
presented in Sect. [2l The initial density was 
expressed via the filling factor (j) = p/po with 
= 0.15 and po = 2000 kg m~^. For the ten- 
sile strength we used the semianalytical relation, 
derived in Sect. | 3.1 | fEq. [S]) t h at ma tches the find- 
ings of iBlum and Schraplerl ( 20041 ) and paper I. 
The bulk modulus was modeled by a power law 



K{p) =Ko^ 



(15) 



with Kq = 300 kPa and the initial density of the 
dust aggregate pi = 300 kg m~^. The bulk modu- 
lus Kq = piC^ for uncompressed material was de- 
termined by the measurement of the sound speed. 



For the compressive strength, several differ- 
ent relations were tested. At first we adopted 
the relation from t h e un iaxial experiments by 
Blum and Schraplerl (|2004l ). Secondly we used 



the omnidirectional compression curve presented 
in this paper. After it turned out that a modified 
relation for the dynamical compressive strength 
curve had to be considered, the omnidirectional 
compression curve (Eq. [10]) was shifted towards 
lower pressure regimes using the free parameter 
Pm (see Sect. 13. 5p . 

Since no experimental data was available for the 
shear strength F, parameter studies were carried 
out with three different rela tions y ( |T|, S ): Y = 
|T|, r = S and, following [Siron^ (|2004f ). Y = 



Y^I]|T|, which represents the geometric mean of 
both quantities. 

Due to reasons of stability, the two materi- 
als in contact (solid projectile, dusty target) have 
to be separated by artificial viscosity. We use 
the approach by onaghan and Gingo 3 (Il983h 
and apply an a- viscosity of 1.0 to all particles of 
the sphere and all particles interacting with the 
sphere. All other dust-sample particles we re sim- 
ulated without artificial viscosity following ISironol 
(|2004j ). All details regarding the SP H code can be 
found in iGeretshauser et al. ( 20091 ). 



4.2. Calibration procedure 

With the aim of reproducing the experimental 
results presented in the laboratory section, an SPH 
simulation using the omnidirectional compressive 
strength curve (ODC) was conducted. In the re- 
sulting pressu re regime the ODC r e lation and the 
relation from um and Schraplerl (|2004l ) are al- 
most identical. Therefore they can be treated as 
one case. 

The impact velocity of the 1.1 mm glass bead 
was 0.65 m s~^ and we will compare the results 
of the simulation with the vertical density pro- 
file along a line through the center of the sphere 
perpendicular to the bottom of the dust sample, 
which was measured with x-ray micro-tomography 
as described in Sect. 13. 4[ Figure [12] shows the 
results for three different shear strength models, 
which are compared to two density profiles as mea- 
sured in the experiments (lines with blue and green 
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crosses). The initial surface of the dust sample is 
at mm. 

For the original ODC relation (pm = 13 kPa), 
the simulations for all shear strength models re- 
sulted in a much too shallow intrusion depth and 
an insufficient maximum filling factor underneath 
the sphere. These findings indicated, that the 
compressive strength curve had to be modified in 
order to reproduce the experimental data. There- 
fore, we performed a parameter study varying 
the parameter i-e. shifting the compressive 
strength curve to lower pressures for the different 
shear strength mod e ls. Fo r the complete study see 
Geretshauser et aP ( 2QQ9I ). Independent experi- 
ments (paper III) also support a lower which 
can quantitatively explain the amount compres- 
sion in bouncing collisions. 

A significant increase of the intrusion depth was 
only observed in case of Y = y^^\T\ and Y = H 
(see Fig. [121 top ^ind center). In case of Y = |T| 
the intrusion depth hardly changed with decreas- 
ing (Fig- US bottom). Since the shear strength 
remained constant and changing did not have a 
significant effect, it can be concluded that shearing 
plays an important role during the intrusion. 

Compared to the other cases, the shear strength 
reaches its highest values in the Y = \T\ case. 
Hence, the material can hardly be pushed away 
due to shear and has to be compressed. Therefore, 
the highest filling factors can be found in this case 
(see Fig. [121 bottom). The Y = H model yields 
the lowest shear strength values. Hence, material 
is mostly sheared aside, less material is compressed 
and therefore this model leads to filling factors be- 
low the reference data (see Fig. [121 center). 

Figure [13] shows intrusion depth over stopping 
time regarding the shear strength models Y = H 
and Y = ^/^T\ for ah p^. Since Y = \T\ did 
only yield insufficient intrusion depths, this model 
was omitted here. A good time/depth match was 
achieved for p^ = 3.9 kPa using F = E and for 
Pm = 1-3 kPa using Y = ^S|T|. However, the 
y = E model cannot reproduce the high values in 
the vertical filling factor profile (Fig. [121 center) 
whereas the Y = y^I;|T| model yields an almost 
perfect match (Fig. [121 top). Therefore, the latter 
with Pm = 1.3 kPa gives a good match in Fig. [12] 
(top) as well as in Fig. [13] and is therefore used 
for further simulations. A more detailed study on 
the determination on the Pm value can be found 




-1.5 -1 -0.5 

intrusion depth [mm] 

Fig. 12. — The filling factor measured along a line 
through the center of the sphere, perpendicular 
to the bottom of the dust sample (lines with blue 
and green crosses). The initial surface of the dust 
sample is situated at mm and the steep slope 
at the right end of the experimental curves marks 
the transition from the dust sample to the glass 
sphere. The other curves are numerical simula- 
tions, varying the shear strength model and the 
material softness p^- The shear strength relation 
was (top) Y = V'EITI, (center) F = S, and 
(bottom) Y = \T\. 
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Fig. 13. — Stopping time - intrusion diagram. Ex- 
perimental data for spheres of 1 mm radius. La- 
bels indicate the pm values (in kPa) for the mod- 
ification of the compressive strength curve. For 
the y = E model the best match in stopping time 
and intrusion depth is found for =3.9 kPa. For 
the Y = Y^I]|T| model the best approximation is 
given for pm = 1.3 kPa. 



in lGeretshauser et al.l ( 2QQ9I ). 

Hereby, we have determined parameters for all 
previously unknown material relations and thus 
have calibrated the SPH model with respect to 
the presented experiments. The resulting strength 
curves of compression (Eq. [TOl Pm = 1-3 kPa), 
tension (Eq. [8]) and shear {Y = y^E|T|) are illus- 
trated in Fig. [Ml 

However, the fact that the filling factor does 
not rapidly drop to ~ 0.15 at a depth of 1.5 mm 
requires further investigation. 

4.3. Reproducing Experimental Features 

Since intrusion time and depth as well as the 
filling factor profile underneath the sphere have 
been used to determine Pm and the correct shear 
strength model, further features have to be repro- 
duced in order to validate the calibration. 

One of these features is the cumulated volume 
over filling factor relation (Fig. [T5|) . While the 
filling factor profile only displays a cut through the 
compressed volume, this curve represents the total 
compressed volume with its filling factors. Both 
curves are not fully, but mostly independent from 
each other. The chosen model and Pm value yield 
an almost perfect match for filling factors > 0.22. 
The deviation for lower filling factors is due to the 
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Fig. 14. — Compressive strength curve S(0) (Eq. 
fTO]) for Pm = 1.3 kPa, tens ile stre ngth |T| (Eq.[8]), 
and shear strength Y = ^/T^^\T\. 
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Fig. 15. — Cumulated volume over filling factor. 
While Pm = 13 kPa, i.e. the unmodified ODC rela- 
tion, and Pm = 0.39 kPa as well as pm = 0.65 kPa 
yield too small and too high compression values, 
respectively, Pm = 1-3 kPa matches very well for 
> 0.22. The experimental data are identical to 
those shown in Fig. [11] (bottom) 



larger amount of compressed volume. This effect 
was already seen in the filling factor profile and 
is also very prominent in the comparison of the 
spatially density distribution plots (compare Figs. 
El and [ED. 

Another feature to be reproduced is the rela- 
tion D (X mvA~^ found in a similar way in the 
drop experiments (cf. Fig. [6]). We performed a se- 
ries of 2D simulations with spheres of 1 mm and 3 
mm diameter and evaluated the maximum intru- 
sion depth with respect to the impact velocity v. 
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Fig. 16. — Spatially distributed compression as 
calculated in the SPH simulation with Y = 
Y^S • |T| and pm = 1-3 kPa; same color scale as 
Fig. [9l the projectile is not plotted. Although the 
filling factor of compressed material is compara- 
ble to the one in the experiments, the compressed 
volume reaches significantly deeper. 
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Fig. 17. — In the momentum-intrusion relation, 
the agreement between simulation and experimen- 
tal results is very good for values of mvA~^ > 
1 kg m~^ s~^. 

The latter was varied from 0.1 m s~^ to 1.0 m s~^ 
in steps of 0.1 m s~^. 

2D simulation and experiment cannot be com- 
pared directly due to the different geometry (the 
2D setup represents a cut through an infinitely 
long cylinder). The advantage of using the quan- 
tity mvA~^ instead of the more accurate Eq. [12] 
is given by the fact that the former can be "con- 




normolized time t' 

Fig. 18. — Normalized deceleration curve com- 
pared to the results. The deceleration curve in 
the SPH simulation is slightly lower than the ex- 
perimentally observed sine curve, but well within 
the errors. This effect will be be analyzed in fu- 
ture work. However, the range of experimental 
data encompasses the simulation results. 

verted" into 2D by the following correction: 

msDV _ ^Trr^p-v _ 8 nr'^p-v _ 8 m2DV 

AsD Trr^ Stt 2r Stt A2D 

(16) 

In comparison with the experimental results, the 
data from the simulation matches very well for 
mvA~^ > 1.0 kg m~^ s~^ (Fig. [T7|). For smaller 
values the simulation yields a shallower intrusion 
than the reference experiments, which, however, 
also show significant scattering in this range. 

Comparing the simulated and experimentally 
acquired normalized deceleration curves (Fig. [18]), 
the simulated data slightly deviate from the exper- 
imental mean but remain within standard deriva- 
tion limits. The deviation could arise from the 
geometric difference of the 2D and 3D case and 
has to be investigated in future works. 

5. Application of SPH to Dust Collisions 
in PPDs, Conclusions, and Outlook 

In this section we will present some prelimi- 
nary applications of SPH simulations to dust col- 
lisions in protoplanetary disks. We will present 
two examples of previously unfeasible calculations 
of inter-particle collisions among macroscopic dust 
aggregates and will qualitatively compare them to 
similar dust experiments performed in the labora- 
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Fig. 19. — Sequence of snapshots of an SPH sim- 
ulation of a fluffy dust aggregate with a radius of 
0.5 mm, impacting a solid target at a velocity of 
0.2 ms~^. The time differences between subse- 
quent images are 0.35 ms, 0.32 ms, 0.23 ms, and 
4.45 ms, respectively. The colors denote different 
degrees of internal compaction. [See the electronic 
edition of the Journal for accompanying mpeg an- 
imations.] 

tory. Then, we will speculate about how the SPH 
code should be used in research on protoplane- 
tary growth. Finally we will sketch future work in 
preparation. 

5.1. Qualitative comparison between SPH 
simulations and laboratory experi- 
ments 

The strength of the SPH simulations - be- 
sides the well-known examples in hyper-velocity 
collisions - over laboratory experiments and 
molecular-dynamics simulations is that low- velocity 
collisions among arbitrary dust aggregates can be 
investigated. Here, we show two examples re- 
cently observed in the lab, which can so far not be 
described by any other model. Example 1 deals 
with the frequently-observed bouncing collisions in 
aggregate- aggregate interactions. Example 2 de- 
scribes the impact of a single dust aggregate onto 
a solid flat target, which shows the co-occurrence 
of (partial) sticking and fragmentation. 



5.1.1. Example 1 

Bouncing in collisions between dust aggregates 
has been observed in many laboratory experiments 
(|Blum and MiinchI Il993l : iLangkowski et all l2008l . 
paper HI; HeiBelmann et al., in prep., will appear 
in this series), although molecular-dynamics simu- 
lations always show a direct transition from stick- 
ing to fragmentation when the collision energy 
exceeds a threshold va l ue ( Dominik and TielensI 
19971 : IWada et al.ll2007l . l2008h . Nature obviouslv 
chooses a wider bouncing transition between those 
two stages, at least for aggregates above a certain 
size. It turns out that the SPH method is capable 
of describing the bouncing phase quite well. We 
have run a 3D SPH simulation of a low-velocity 
impact of a 1 mm (diameter) fluffy aggregate onto 
a flat target. Due to symmetry arguments, this 
is identical to a two-aggregate (central) collision 
with twice the collision velocity. In our case, the 
aggregate was composed of 33,377 SPH particles 
and had an initial volume filling factor of 0.15. 
All other material parameters were identical to 
those in the previous section, i.e. Kq = 300 kPa, 
= 1.3 kPa, Y = ^/^T\. The impact velocity 
was 0.2 ms~^, matching exactly the situation in 
the aggregate-wall experiments performed in pa- 
per HI and also those in the aggregate- aggregate 
collisions investigated by HeiBelmann et al. (in 
prep.) with a collision speed of 0.4 ms~^. Fig. [19] 
shows a sequence of snapshots with a cut through 
the center of the aggregate, indicating the inter- 
nal compaction due to the impact. Our simulation 
can corre ctly predict the coefficie nt of restitution 
of - 0.2 (|Blum and MiinchI [19931 . HeiBelmann et 
al., in prep.), although details in the compaction 
behavior still deviate from the laboratory results, 
which might be caused by insufficient resolution in 
the SPH simulation. 

5.1.2. Example 2 

In the previous example, we have seen that 
bouncing marks the broad transition regime be- 
tween sticking and fragmentation. However, in 
the case of the impact of a dust aggregate onto 
a solid target, laboratory experiments have shown 
that, for impact experiments above the fragmen- 
tation threshold, fragmentation is always accom- 
panied by partial sticking of the agg regate to the 
target. This effect was first found bv I Wurm et al ' 
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Fig. 20. — Image sequence of an experiment, in 
which a fluffy dust aggregate impacts a sohd tar- 
get at 8.4 m s~^. Part of the aggregate sticks to 
the target after the cohision. [See the electronic 
edition of the Journal for accompanying mpeg an- 
imations.] 




Fig. 21. — Image sequence of an SPH simula- 
tion, identical to the experiment shown in Fig. [20l 
The simulation with the calibrated parameters 
(left) cannot reproduce the experimental results, 
whereas a simulation using an unidirectional com- 
pression curve and E = |T| (right) can repro- 
duce the qualitative findings of Fig. [20l [See the 
electronic edition of the Journal for accompanying 
mpeg animations.] 



(|2QQ5l ) for compacted dust aggregates and impact 
velocities above 25 ms~^ and later confirmed in 
our laboratory for <j) = 0.35 aggregates and im- 
pact velocities above 1 ms~^. Fig. [2Ql shows an 
image sequence of an impact experiment with frag- 
mentation and partial sticking. An average of 
10 % of the projectile mass sticks to an initially 
smooth target at normal impact, wh ich is consis- 
tent w ith the low velocity results of IWurm et al.l 
([20051). The remainder of the projectile mass 
is fr agmented into a power- law mass distribution 
(see Blum and Miinchll 19931 ). The fragments leave 
the target under extremely flat angles. Our SPH 
simulation (Fig. [211 l^ft) featuring the calibra- 
tion parameters of Sect. 14.21 cannot reproduce the 
fragmentation behavior seen in the experiments. 
Here, the predominant part of the dust sample 
sticks to the target. Only a few bigger chunks 
and single SPH particles burst off. However, a 
simulation with the same setup, but using the 
shifted unidirectional compressive strength rela- 
tion (Sect. 13. 2p and a shear strength that is equal 
to the tensile strength, matches the experimental 
observations at least qualitatively (Fig. [211 right). 
From that we conclude that the SPH code is in 
principle capable of simulating fragmentation of 
highly-porous aggregates, even without the dam- 



age model adopted in the original iSiron 3 (|2004l ) 
porosity model. 

We conclude that the shear model Y = S^-^ • 
|T|^-^ tested for the dynamic compression exper- 
iments (Sect. 14. ip is unable to explain the frag- 
mentation findings which are rather dominated by 
shear and tension, whereas a shear model Y = 

• |T|^ shows qualitative agreement. The imper- 
fect shear model can also be responsible for the 
narrow but deep compressed volume in Fig. [16] 
compared to Fig. [9l A future task will therefore 
be to refine the shear calibration in a way that we 
will use y = C • • |T|^~" with the free param- 
eters C and a. Comprising both experiments for 
calibration we will be able to find a shear model 
that can reproduce both cases. 

5.2. Use of the SPH code in research on 
protoplanetary growth 

The above examples show that the SPH method 
is a powerful tool to investigate the outcomes of 
protoplanetary dust collisions. When properly cal- 
ibrated with laboratory experiments, SPH calcu- 
lations allow access to parameter-space regions 
that are unavailable to laboratory experiments. 
Whereas molecular-dynamics simulations can be 
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used for studying collisions of very small dust ag- 
gregates, SPH is most useful for very large sam- 
ples. Such samples, particularly those with fluffy 
compositions, cannot be built or treated in labo- 
ratories, and the experimental study of collisions 
seems impossible. 

A particularly interesting and still unsolved 
problem is the dichotomy in the collision behavior 
of pairs of dust aggregates with similar and differ- 
ent sizes, respectively. In paper II we found stick- 
ing by deep penetration for impacts of mm-sized 
dusty projectiles into flat, cm-sized dusty targets 
("projectile-target" collisions) above ~ lms~^. 
Both dust aggregates, projectile and target, con- 
sisted of identical particles and had equal poros- 
ity. Using similar dust aggregates, but giving 
projectile and target comparable siz e ( "pro j ectile - 
projectile" collisions), Blum and Miinch ( 19931 ) 
and HeiBelmann et al. (in prep.) found that colli- 
sions either lead to bouncing or to fragmentation. 
Bouncing instead of sticking was also observed 
in paper II when the target aggregates were pre- 
pared such that the local radius of curvature cor- 
responded to the projectile's radius. To find out 
where the boundary between "projectile-target" 
and "projectile-projectile" collisions occurs, will 
be one of our future applications of our SPH code. 

5.3. Future work 

We have only begun to explore the potentials of 
SPH simulations of collisions between protoplane- 
tary dust aggregates. Before we can start to inves- 
tigate the full parameter space in protoplanetary 
dust collisions, i.e. before we can begin to find out 
what the collisional outcome is for all combina- 
tions of aggregate size, porosity, collision velocity, 
impact angle, state of rotation, temperature and 
state of sintering, material and size (distribution) 
of the constituent dust grains, etc., the material 
parameters of macroscopic dust aggregates have 
to be fully explored. This will be the next task in 
our investigation. To achieve this, we will perform 
more calibration experiments of the type described 
in this paper for dust aggregates of various com- 
positions and porosities. In addition to that, other 
calibration experiments will be explored, like the 
ones described in Sects. 15.1.11 and 15.1.21 
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